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Abstract 

Considering of the possible difference in strange and antistrange quark 
distributions inside nucleon, we investigate the — Dj asymmetry in pho- 
toproduction in the framework of heavy-quark recombination mechanism. We 
adopt two distribution models of strange sea, those are the light-cone meson- 
baryon fluctuation model and the effective chiral quark model. Our results 
show that the asymmetry induced by the strange quark distributions is dis- 
tinct, which is measurable in experiments. And, there are evident differences 
between the predictions of our calculation and previous estimation. Therefore, 
the experimental measurements on the — D~ asymmetry may impose a 
unique restriction on the strange-antistrange distribution asymmetry models. 
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High energy charmed hadron production plays an important role in studying strong 
interactions. Due to the large charm quark mass, the charm quark involving processes 
can often be factorized, i.e., into a hard process, which describes the interaction scale 
much larger than Aqcd, the typical QCD renormalization scale; and a soft part, which 
demonstrates the hadronization effects. In recent years, a large production asymmetries 
between the charmed and anti-charmed mesons have been measured in fixed-target hadro- 
and photo-production experiments [TJ El HI IB El EI • Among these, photoproduction is 
thought to give more clean signature than the hadroproduction ones, because there is only 
one hadron in the initial state. However, it is intriguing to notice that the experimental 
observation on the asymmetries of the charmed hadron production are greatly in excess 
of the predictions of perturbative QCD(pQCD) [HI GDI HQ • 

For large transverse momentum processes, the QCD factorization theorem en- 
ables us to calculate the cross sections of heavy hadron production by the pQCD. To be 
specific, for the D s photoproduction, the differential cross section may be expressed as 
the convolution of the parton distribution function, the partonic cross section, and the 
fragmentation function, like 

d<7 [ 7 + N^D a + X\ = Y,fi/N® Ml + c + c + X] (g) D^ D „ (1) 

where /j/jvis the distribution function of the parton % in a nucleon iV; da (7 + 1 — > c+c+X) 
is the pQCD calculable cross section of partonic subprocess; and D c ^c s is the nonper- 
turbative fragmentation function. In this picture, c and c are produced symmetrically at 
leading order in a s . The asymmetry appears only in the next-to-leading order (NLO), 
or higher, corrections. However, the charm- ant icharm asymmetries predicted by NLO 
correction are an order of magnitude smaller than the asymmetries observed in photopro- 
duction experiments |13j . That is, the charmed and anticharmed photoproduction asym- 
metry may not be fully explained by the charm- ant icharm quark production asymmetry. 
Moreover, the theoretical predictions through above mechanism (J!Q) cannot account for 
the differences among the asymmetries of D mesons with different light quark flavors. 

There have already been some early attempts to explain the observed asymmetries 
Ho] . In these approaches, the asymmetry is supposed to appear due to the nonperturbative 
hadronization effects. Hence, they are all sensitive to unknown distribution functions of 
partons in the remnant of the target nucleon or photon after the hard scattering. In 
comparison, the heavy-quark recombination mechanism proposed by Braaten et al. can 
give a more reasonable explanation to the D meson and A c production asymmetries. In 
the heavy-quark recombination mechanism, a light parton(g) that participates in the hard- 
scattering process recombines with a heavy quark(c) or an antiquark(c) and subsequently 
hadronize into the final-state heavy-light meson. The recombination happens only when 
the light-quark in the final state has momentum of O(Aqcd) in the heavy-quark, or 
antiquark, rest frame. Namely, the light-quark and the heavy-quark recombine in a small 
phase space. By virtue of the heavy quark symmetry, SU(3) flavor symmetry and the 
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large N c limit of QCD, the heavy-quark recombination mechanism gives a simple and 
predicative explanation for the asymmtries with two nonperturbative parameters. Since 
the light quark q in the recombination model can be either u, d or s quark, it can account 
for the difference of asymmetries among different light flavors. 

In Ref. [Ej, Braaten et al. calculated the charm/anticharm production ratio and 
asymmetry of D mesons using the heavy-quark recombination mechanism and confronted 
their results to the experimental data. In their consideration, the asymmetry of D s meson 
comes from the process in which the c(c) and light valence-quark of nucleon recombine 
into a D~(D + ) meson, while the recoiling c(c) quark fragments to Df(D~) meson. That 
is, the — D~ asymmetry stems from the excess of u and d over u and d in the 
nucleon. Because the s and s content of the nucleon are identical in their assumption, 
the asymmetry of D s meson has the opposite sign as that of D + and D~ meson and is 
relatively small. The experimental data on D s exist [3|7j, but with very large errors, and 
hence do not tell whether there is an asymmetry or not. Different from their consideration, 
recent years, many studies show that there is striking strange/antistrange sea asymmetry 
in the momentum distribution inside the nucleon [T%1 H§] . Stimulated by this idea, in 
this work, we calculate the D s meson production asymmetry induced by the strange- 
antistrange quark distribution asymmetry within nucleon by employing the heavy-quark 
recombination mechanism. 

According to the heavy-quark recombination mechanism, the cross section of D s meson 
photoproduction may schematically formulated as: 

da [y + N^D S + X] = f q/N ®Y1 <®h + * - M n + c]p[(cs) n - D s ], (2) 

where (cs) n represents that the s in the final state has small relative momentum in the c 
rest frame, and n is the color and angular momentum quantum numbers of (cs) interme- 
diate state. da[j + s — > (cs) n + c] is the perturbative QCD calculable partonic subprocess. 
The factor p[(cs) n — > D s ] is the probability of the (cs) n state to evolve into a final state, 
here, the D s . 

In our consideration the D s meson may be produced via two different schemes, i.e. 

(a) daft + (cs) n + c}p[(cs) n -> D s ] , (3) 
(ft) + q - (cq) n + c}J2 p[fa) n - D] ® D„ D . , (4) 

D 

while in Ref. |17j only the second one, the recombination-fragmentation mechanism, was 
taken into account. In process (a), the (cs) n recombines into the D s meson directly; in 
process (b), (cq) n recombines into the Dj, D~ or D° meson, and the recoiling c quark 
fragments to the meson. 

The calculation of the partonic cross section da[y+q — > (cq) n + c\ in pQCD is straight- 
forward, and our calculation confirms the Ref. ^7j results. That is, 

d& i- Ac(i)\i 2567r 2 e 2 c aa 2 s m\ r S kT 2 

~di [cq[b ° )l ~ 81 ^ h f7 (1 + T } 
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mlS, S 3 2(1 + k)S a k 2 T^ 2mlS 3 kT 

+ T + + + (5) 



m^S 3 4(2 + k)S 2 2(3 + 7k)^ , tn , 3k 2 T, 

+ -^(y3+ T2 + r +4«(3 + «) + — ) 

Here, k = e q /e c is the ratio of the electric charge fractions of light and charm quarks. The 

2 T 1 — I rn m \2_^2 „„J TT — _^ \"i_^2 
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Lorentz invariants are defined as S = (p q +p-y) 2 , T = (p 7 — p c ) 2 — m 2 , and U = (p 7 — p. 
p q , p 7 and p c are the momenta of the light quark, photon, and c quark, respectively. It 
is noted that because the relatively small momentum of the light quark in (cs) n system, 
the higher angular momentum excited states are suppressed by power of ^qcd/pt or 
KqcdI^c- For a leading order estimation, we take only the contributions from 1 Sq and 
3 Si states. 

To calculate the total cross section (da[j + N — > D s + X]) of the production of D s 
meson, we need to know the distribution of strange sea quark. According to the light- 
cone meson-baryon fluctuation mo del [IB], the strange quark-antiquark asymmetry of the 
nucleon sea is generated by the intrinsic strangeness fluctuations in the proton wavefunc- 
tion. In this model, the asymmetry stems from the intermediate K + A configuration of 
the incident nucleon, which has the lowest off-shell light-cone energy and invariant mass 
[20j . The momentum distributions of the intrinsic strange and antistrange quarks in the 
K + A system can be parameterized as 

8 (x) = I —f\/K+\{y)q 3 /A{-) , (7) 

J x V V 



s( x ) = [ —f K +/ K +A(y)qs/K+(-) , (8) 
J* v y 

where f k/ k+ k(jj\ fK+/K+\(y) ar e the probabilities of finding A, K + in the (K + A) state; 
9s/a(|) ; Qs/k+(^) are probabilities of finding strange and antistrange quark in A or K + 
states. To estimate these quantities, two simple functions of the invariant mass M 2 = 
S 2 = i K±i ^ mi f° r the two-body wavefunction are given [2*Tj . 

4>Gaussian(M 2 ) = A Gaussian eXp{-M 2 /2ot 2 ) , (9) 

^ PoW er(M 2 ) = A power (l + M 2 /a 2 )~ p , (10) 
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where the a represents the typical internal momentum scale. In our analysis in this work, 
we simply adopt the Gaussian type wavefunction for use. 

In recently, another model, which based on the effective chiral quark theory, for the 
distribution of strange sea quark is proposed by Ma et al. [I9j. In this model, the strange 
quark distribution is determined by both the constituent quark distribution and the quark 
splitting function. For instance, 

S{x) = P s K+/u<S> U + PsK0/d<S>d , (11) 

and 

8(X) = V ~s/K+ ® P K+s/u (g) U + V s/K° (g) P K0 S /d (g) t*0 • (12) 

Here, Pj as /i is the splitting function representing the probability of finding a constituent 
quark j together with a spectator Goldstone(GS) boson (a = tt,K,t]) in a parent con- 
stituent quark. w and do denote the constituent quark distributions; Vi/ a is the quark i 
distribution within the GS boson a. 

Using the resultant production cross section, it is straightforward to calculate the 
asymmetry of the — D~ photoproduction, which are defined as 

a[D s ] = °* ~ a °° . (13) 

In our calculation, the two nonperturbative input parameters, p sm and p s f, are extracted 
from experiments by fiting to the E687 and E691 data [S1IZ1- 

Psm = Peff[cu( 1 S )^D }=p eff [cd( 1 So)^D + }=p eff [cs( 1 S )^D+} 

= Peff [cu( 3 S l ) - D*°] = Peff [cdCS l ) - D*+] = p eff [cs{ 3 S 1 ) - D*+) = 0.15 , 



(14) 



and 



Psf = Peff[cu( i S 1 ) - D°] = PefficdCS,) - D + ] = p ef f[cs( 3 S 1 ) - D+] 

= Peff [cuCS ) D*°] = Peff [cdCS ) - D*+] = p eff [cs( l S ) ^ D*+] = . 

(15) 

Here, the subscripts sm and sf stand for spin-matched and spin-flipped situations, re- 
spectively. Considering the possible uncertainties, the above values are in agreement with 
those extracted from the hadroproduction 



As for the fragmentation function of charm quark to the D s meson, we exploit the 
well-known Peterson fragmentation function, 

D c ^ Ds (z)=P(z;e)f c ^ Ds . (16) 
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Figure 1: The inclusive cross sections da/dxp, calculated by using the light-cone meson- 
baryon fluctuation model, for D s production with p sm = 0.15 and p s f = 0. The solid, 
dash-dotted, and dotted lines correspond to the production of Df, D~ in recombination 
mechanism, and Df at the leading order photon-gluon interaction with the fragmentation 
mechanism, respectively. 



Here, D(z; e) j2H| is the Peterson function and / c _» n is the fragmentation probability. 
In practice, the parameter e is set to be 0.06, as in Ref. ^7j, which falls in the region of 
experimental fitting allowance. All through this paper, the charm quark mass are set to 
be 1.5 GeV, and the factorization and renormalization scales are set to be at the meson 
transverse mass mx, the \Jp\ + mpj s - The parton distribution distribution function(PDF) 
of CTEQ6L 25j is used. For simplicity, we use the average photon beam energy of the 
E687 experiment, (E^) = 200GeV, in the calculation. To be more accurate, the D* feed 
down effect is included in our analysis, but those from other higher excited states are 
thought to be small and simply neglected. 

In figures 1 and 2, based on the light-cone meson-baryon fluctuation model, the depen- 
dence of the inclusive cross sections of Df and D~ production on the Feynman variable 
Xp and the rapidity y respectively are shown. The corresponding effective chiral quark 
model results are presented in Figures 3 and 4. For comparison, the cross section of the 
leading order photon-gluon fusion process, which will give a symmetric Df production 
through fragmentation, is also presented. We find that the results appear to be insensitive 
to the variations of the factorization scale and the Peterson parameter e. From Figs. 1 - 
4, we see that cross sections peak in the region of small xf and large rapidity y. Although 
the contributions from the recombination mechanism are not the dominant ones, they 
may give out the asymmetry. Both for the light-cone meson-baryon fluctuation model 
and the effective chiral quark model, in the region of the experimental cut xp > 0, the 
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Figure 2: Inclusive cross sections da/dy calculated with the light-cone meson-baryon 
fluctuation model. The curves are described as in Figure I. 
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Figure 4: Inclusive cross sections da/dx y calculated with the effective chiral quark model. 
The curves are described as in Figure 1. 

corresponding cross sections for are bigger than those for the D~. Thus, the D s pho- 
toproduction asymmetry appears with the effect of strange-antistrange quark asymmetric 
distributions. 

Our predicted D s production asymmetries, in comparison with Ref.|17j. are shown in 
in Figure 5. We find that, by adopting the light-cone meson-baryon fluctuation model, 
the production asymmetry of D s meson is about 1.2 times larger than that in Ref . 11 7] ; 
and by using the effective chiral quark model, it is about 80% bigger only. If we extend 
to the xp negative region, the asymmetries coming from the intrinsic strange sea momen- 
tum distribution are very big and diverge very much from the prediction of Ref.^Tj. In 
this region, the asymmetry from the light-cone meson-baryon fluctuation model is even 
flipped relative to the positive Xp region, while the effective chiral quark model and the 
fragmentation-recombination scheme give the results with the same sign as in the Xp > 
region. Since the obtained data are in the xp > region, here we will not show the results 
in the xp < region. From the figure 5, it is obvious that the three different asymmetry 
producing schemes(models) give results diverging with each other with the xp increase. 
This scaling difference leaves the experiments with an opportunity to decide the physical 
reality. 

In summary, in this paper, the heavy-quark recombination mechanism was employed, 
which gets a first success in explaining the D meson and A c asymmetry production. We 
have studied the — Dj asymmetry in the photoproduction. Our point is that this 
asymmetry can be induced not only by the recombination-fragmentation mechanism, the 
Eq01 but also by the strange and antistrange distribution asymmetry inside the nucleon, 
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Figure 5: The asymmetry a[D s ] versus xp. The dotted and dash-dotted lines correspond 
to the results from the light-cone meson-baryon fluctuation model and the effective chiral 
quark model, respectively. The solid line from the Ref.[T7] result. 



the EqO And, we find the latter effect is even bigger then the former one depending on the 
employed model and phase space region. Two QCD relevant strange quark distribution 
models were used, that is the light-cone meson-baryon fluctuation model and the effective 
chiral quark model. After including the process 0, the predicted asymmetry increases 
anyhow. However, what we discussed should be the dominant contributions to the D s 
photoproduction asymmetry. Although there remains some uncertainties, such as from 
the breakdown of SU(3) symmetry, the use of large N c limit and the NLO corrections to 
the leading order photon-gluon fusion, etc., our results are adequate to make qualitative 
predictions, and more detailed discussion on the relevant uncertainties can be found in 
Ref.|lfil ITT]. It is noticed that the nowadays D s asymmetry experimental data are very 
limited and preliminary. With enough data collection in the future, it is expected that 
the experiment can inversely impose a strong restriction on the strange and antistrange 
quark distributions by measuring the D s production asymmetry. 
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